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ABSTRACT 

We use high resolution N-Body/SPH simulations to study the hydrodynamical and 
gravitational interaction between the Large Magellanic Cloud and the Milky Way. We 
model the dark and hot extended halo components as well as the stellar/gaseous disks 
of the two galaxies. Both galaxies are embedded in extended cuspy LCDM dark matter 
halos. We follow the previous four Gyrs of the LMC's orbit such that it ends up with 
the correct location and orientation on the sky. Tidal forces elongate the LMC's disk, 
forcing a bar and creating a strong warp and diffuse stellar halo, although very few 
stars become unbound. The stellar halo may account for some of the microlensing 
events. Ram-pressure from a low density ionised halo is then sufficient to remove 
1.4 X IO^Mq of gas from the LMC's disk forming a great circle trailing stream around 
the Galaxy. The column density of stripped gas falls by two orders of magnitude 100 
degrees from LMC. The LMC does not induce any response in the Milky Way disk. 
On the contrary, the tides raised by the Milky Way determine the truncation of the 
satellite at about 11 kpc. After several Gyrs of interaction the gas disk of the LMC is 
smaller than the stellar disk due to ram pressure and its size compares well with the 
observational data. 
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1 INTRODUCTION 

The Large Magellanic Cloud (LMC) is in the process 
of destroying its neighbour, the Small Magellanic Cloud 
(SMC). In turn, the Milky Way (MW) is playing a role 
reshaping these galaxies and creating the spectacular Mag- 
ellanic Stream (hereafter MS). The MS is a trailing filament 
of neutral hydrogen that originates from the Magellanic 
Clouds and stretches for over ~ 100 degrees in the Southern 
Sky and is the most prominent signature of an interaction 
between the Clouds and the Milky Way. The combined 
interaction between these three galaxies must strongly 
influence the morphological evolution of the Clouds, as well 
as their internal kinematics and star formation history. 
A resonant interaction between the Galaxy and the LMC 
could even be sufficient to excite th e warp observed in the 
Galactic disk ^Weinberelll995l ■ ll998^ . 

Several models have been proposed in order to study the 
nature of this interaction and to isolate the main mechanism 
responsible for the formation of the Stream. These models 
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have become increasingly refined over the years. The first 
gravitational tidal models took into acc ount Just the binary 
inter a ction between the LMC a nd MW iLin fc Lvnden-BelJ 
llQTTt iMurai fc Fuiimotolll980H . but later models included 
a close interaction between the Clouds themselves during a 
previous passage at the perigalacti c on about 1.5 Gyrs ago 
'Gardiner. Sawa, & Fuiimoto 199J; |l in. Jones, fc Klemolal 
1995; Gardiner & Noguchi 1996). 

The strongest argument which supports a tidal origin 
for the MS is provided by the discovery of n eutral hydro- 
gen in front of the LMC jPutman et aljll998h . However the 
amount of mater ial in this feature is j ust a fraction of the 
mass of the MS JPutman et alj|2003all . However there are 
some characteristics of the Stream that tidal models cannot 
explain. In particular they fail in reproducing the gradual 
decrease of the HI column density and the lack of a stel- 
lar tidal feature ijMathewso^^ aL^979 HBrueck fc Hawkins! 
119831: iGuha thakurt a fc Reitzeill99ar Stars are expected to 
be strip ped by gr a vitatio nal forces along with the gas, al- 
though IWeinbera i200(il suggests that the stellar debris 
would have a different distribution respect to the gaseous 
Stream and form a diffuse envelope around the LMC. Traces 
of an offset old steUar stream could be identified with the 
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giant stars observed by iMaiewski et alJ il999ft around the 
Magellanic Clouds. 

A hydrodynamical interaction has a l so bee n invoked 
to explain the Stream. iMathewson et ahl (^S^) proposed 
that gas is swept away from the intra-cloud region and de- 
celerated by random collisions with high velocity clouds 
(HVC). This model requires unrealistic densities for the 
HVCs in the MW halo in order to justify the high num- 
ber of encounters that should generate the Stream and does 
not explain t he column density gradient. Cont i nuous strip- 
ping models i Mcurcr, BickneU fc Gingoldlll985l:rSofuelll994l: 
iMoore &: Davis J.994) propose a scenario in which gas is 
stripped by ram-pressure forces from the outer regions of 
the Clouds through an interaction with a distribution of hot 
gas in the halo or an extended ionised disk. A satellite mov- 
ing with a relative velocity v respect to this external medium 
of density p experiences a pressure P = pv'^ which can be 
sufficient to remove part of the neutral gas from the disk. 

Together with the difficulties in creating a continuous 
leading arm - although the stripped material decelerated 
by ram pressure forces could fall on smaller orb its and 
lead the Clouds JSofuelll994l: iMoore fc Da^ll994l) - ram 
pressure models are thought to be able to reproduce the 
large amount of gas observed in the Stream (~ 2 x 10^ Mq 
according to Putman et al. 2003a). 

The presence of an extended hot halo surrounding the 
MW and in hydrostatic equilibrium within the dark matter 
potential is expecte d by current models of hierarchical 
structure formation jWhite fc FrenkI Il991^ . This gaseous 
halo, left over from the initial collapse, would be contin- 
uously fuelled by accretion of satellite galaxies over time. 
Recent observational results have confirmed the existence 
of a distribution of hot (T ^ 10® K) gas well beyond the 
Galactic disk, with a radius larger than 70 kpc in order 
to explain some ioni s ation features discover ed in the MS 
JSembach et all 120031: iPutman etaP |20033). Constraints 
from dynamical and thermal arguments fix the density of 
the gaseous halo in a range between 10~^ and lO^^cm"^ at 
the LMC distance from the Galactic center. The Clouds are 
therefore subjected along their orbits to the ram pressure 
generated by a tenuous distribution of hot gas and the 
Magellanic System itself seems to be the result of a complex 
interaction that involves both hydrodynamical and tidal 
processes. 

The aim of this work is to study for the first time the 
simultaneous effect of gravitational and hydrodynamical 
forces acting on the LMC as it moves within the Galactic 
halo, using fully self consistent galaxy models. In particular 
we are interested in the formation and evolution of the 
Magellanic Stream and in the dynamical changes in the 
internal structure of the LMC due to the interaction with 
the MW. We ignore the SMC owing to its small mass 
(~ 10% of the LMC), neglecting the possibility of a close 
encounter between the two Clouds and the consequence of 
additional loss of mass. This approximation would lead to 
an underestimation of the amount of HI in the Stream. 

This paper is structured as follows. In Section 2 we 
present the main characteristics of the galaxy models that 
we use. The results of test wind tube simulations are dis- 



Table 1. Galaxy models: for each component mass (in units 
of IO^^Mq and scale radius (kpc) are indicated. The last three 
columns refer to high resolution runs and report number of par- 
ticles, mass particle (in units of IO^Mq) and softening length 
(kpc). 
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Gaseous disk 
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DM halo 
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Hot halo 
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Mass 
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Gaseous disk 


0.11 


1.66 




10^ 


0.11 


0.1 


Stellar disk 


0.11 


1.66 
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X 10* 


0.22 


0.1 


DM halo 


2.38 


6.02 


6 


X 10^ 


10.4 


0.5 



cussed in Section 3. In Section 4 we analyse the results of 
the interacting runs regarding the Stream and the LMC disk. 



2 GALAXY MODELS 

The initial conditions of the sim ulations ar e cons tructed 
using the technique described bv iHernauisti il993ll . Both 
the Milky Way and the Large Magellanic Cloud are multi- 
component systems with a stellar and gaseous disk embed- 
ded in a s pherical dark matter halo. The density profile of 
the NFW jNavarro. Frenk fc Whitelll997l) halo is adiabati- 
cally contracted due to baryonic cooling dSpringel fc Whit3 
[l99i). Stars and cold (T = 10*K) gas in the disk follow the 
same exponential surface density profile of the form 
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where Aid and Rd are the disk mass and radial scale length 
(in cylindrical coordinates) respectively, while the thin ver- 
tical structure is characterised by the scale height Zd ~ ^Rd- 
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The MW model comprises of two further components: a stel- 
lar bulge and an exte nded hot gaseous halo. The bulge has 
a spherical Hernquist (lHernauislill990l) profile: 



Pb{r) 
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where Mb is the bulge mass and rb its scale length. 

Although the presence of an extended corona of hot 
gas seems to be the most likely explanation for the de- 
tection of OVI and OVII absorption lines associated with 
HI structures (MS, some HVCs, Outer Spiral Arm, Com- 
plex A, Complex C, according to Sembach et al. 2003), 
its density structure is still uncertain, particularly at large 
distances from the Galactic center. These ionization fea- 
tures suggest a mean density of the gaseus halo within 80 
kpc o f ~ 10"* - lO'^'cm"^ (■Sembach.200a' l. .Stanimirovl3 
(|2002) provides similar values for the density at 45 kpc 
{nh ;$3 X 10~*cm~^) in order to explain the confinement 
mechanism of the MS clouds through pressure support from 
the hot halo. Another constraint comes from the detection of 
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Ha emission along the MS which cannot be e xplained solely 
by ph otoionisation from the Galactic disk jPutman et alJ 
I2OO31J1 and could be induced by the interaction with a rel- 
ative dense (uh > lO'^^cm"^ at 50 kpc) gaseous environ- 
ment JWeiner fc WiUiamJl99^ With a density higher than 
10~'^cnr^^ TQinh^^Moora " (l200lh are able to reproduce, us- 
ing hydrodynamical simulations, the head-tail structure and 
the compression fronts associated with HVCs, in the case of 
either pure gas clouds or gas embedded within a dark matter 
halo. A significantly lower value for the upp er limit of the 
halo density (n^ < 10"^cm"^) is found by iMuralil ll2000ll 
by requiring the survival of the cloud MS IV for 500 Myrs 
against evaporation. However we emphasise the fact that 
the evaporation tim e strongly depends on the temperature 
of the gaseous halo (|Cowie fc McKeelll977^ and a choice of 
T ~ 10® K, three times lower than the value used by Mu- 
rali, increases by almost a factor fifteen the lifetime of MS 
IV, making the survival of the cloud possible at a density 
of even 10~*cm~^. At larger distances from the Galactic 
center constraints on th e halo density are poor, although 
iBlitz fc Robishawl lj200nh derive a value of 2.5 x 10"^cm"^ 
at a distance of 200 kpc by applying ram pressure stripping 
to dwarf spheroidal galaxies within the Galactic halo. 

We model the hot halo with a spherical distribution of 
gas which traces the dark matter profile within the virial ra- 
dius and is in hydrostatical equilibrium inside the Galactic 
potential. Assuming an isotropic model, the halo temper- 
ature at a given radius r is determined by the cumulative 
mass distribution M{r) of the dark, stellar and gaseous com- 
ponents of the MW beyond r and by the density profile ph{r) 
of the hot gas: 



300 



T(r) = 



1 



ks Ph(r) 



Ph{r) 5^dr, 
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where mp is the proton mass, G and fcs are the gravitational 
and Boltzmann constants. The halo gas is completely ionised 
at a mean temperature of 10® K. We normalise the density 
profile using a hot halo mass of Mh = 10^° M©, which cor- 
responds to a mean density of 2 x 10~^cm~'' within 150 kpc 
and 8.5 x 10"^ at 50 kpc. With the assumption that the gas 
traces the dark matter we overestimate by a factor of three 
the ionised gas density near the disk, however this is not 
important for the current study. 

The masses and the scale lengths of the MW and 
LMC models (Table 1) are select ed in orde r to r epro- 
duce observational constraints from iKim et aH lll998h and 
Ivan d or Marcl ot al. (2002), wh ile for the MW we adopt a 
model similar to model Ai of iKlvpin. Zhao fc Somervilld 
i2002h . Fig. 1 and Fig. 2 show the rotation curves of the 
different components of the two galaxies. The concentra - 
tion c and the spin parameter A llMo. Mao fc Whitelll99dl 
are indicated at the top left. The concentration is defined 
as cc — r^ir/rs, where r^ir and are the virial and the 
scale radius of the NFW halo, whereas the spin parame- 
ter relates the angular momentum J and the total energy 
E of an halo with a virial mass M^ir through the rela- 
tion A = J\E\^^^G~'^M~^J^. The disk of the satellite galaxy 
has, within the scale radius, a central mass surface density 
of ~ 70Mqpc~'^, which corresponds to a B-band surface 
brightness of ~ 23 mag arcsec"^, if we adopt a mass to light 
ratio of two. 
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Figure 1. MW rotation curves. On the top left spin and concen- 
tration parameter are indicated. 
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Figure 2. LMC rotation curves. On the top left spin and con- 
centration parameter are indicated. 



The thick ness of the disks is set such that Toomre's 
jToomrelll964) stability criterion is satisfied, which requires, 
for a stellar disk 



n ( ^ — ^ 1 



(5) 



where gr is the radial velocity dispersion, k is the local 
epicyclic frequency and Eg the stellar surface density. For a 
gaseous disk Equation [^becomes: 



Q gas (^) 
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where Va is the sound velocity and Eg the surface density of 
the gas. 

AU the simulations we now discuss were carried out 
using GASOLINE, a paral lel tree-code with multi-stepping 
JWadslev. Stadel fc Quinnl l2004l . The high resolution runs 
have 2.46 x 10® particles, of which 3.5 x 10^ are used for 
the disks and 5 x 10^ for the hot halo of the MW (Table 
1). The gravitational spline softening is set equal to 0.5 kpc 
for the dark and gaseous halos, and to 0.1 kpc for stars and 
gas in the disk and bulge components. 



3 TEST SIMULATIONS 

The mass resolution of the simulations strongly influences 
the ram pressure process. In particular, if the particles of 
the hot external medium are too massive, they do not flow 
smoothly on to the disk, but produce scattering and numer- 
ical holes, increasing the stripping process which is domi- 
nated by noise. In order to study these effects, we perform 
several test simulations in which we place the model LMC 
galaxy within a "wind tunnel" and vary the mass of the hot 
halo particles for a fixed resolution of the LMC disk. 

The first set (HI, H2 and H3) of runs simulates the 
passage of the satellite at the perigalacticon, which is the 
point along the orbit where ram pressure is most efficient due 
to the high values of pg and v. The volume of the simulations 
is an oblong of base equal to the diameter of the LMC and 
height h = vt, where v is the velocity of the satellite at 
the perigalacticon (250kms~^) and t is set ~ 2 Gyrs (~ the 
inferred orbital time for the LMC) . We represent the hot gas 
as a flux of particles moving with a velocity v along the z 
axis, parallel to the angular momentum of the disk. The hot 
particles have an initial random distribution, a temperature 
r = lO^K and a number density uh = 8 x 10"^cm"^. The 
box has periodic boundary conditions in order to restore the 
flow of hot gas that leaves the oblong. The disk particles have 
the same mass mdisk as in the high resolution interacting 
runs (see Table 2), while the mass ratio with the particles 
in the oblong mh/mdisk of 1:10, 1:2 and 1:1. Gas is removed 
from the disk if the ram pressure produced by the hot halo 
of the MW is greater than the restoring force per unit area 
provided by the disk of the satellite. The condition for ram 
pressure stripping is expressed by ICunn fc Gott.il972M'l 

Phv^ > 27rGE(ii)E9(i?), (7) 

where i; is the velocity of the galaxy with respect to the sur- 
rounding medium, ph is the density of the hot halo of the 
MW and Eg(r) is the cold gas surface density at the radius 
R. E represents the gravitational surface mass density of the 
disk. The minimum radius given by Equation (7) is the flnal 
stripping radius Rstr beyond which the ISM can be removed 
from the galaxy. The halo component does not contribute 
to the restoring force initially, since in the case of a face 
on impact any force associated to a spherical potential is 
perpendicular to the wind pressure. As the gas is stripped 
out from the disk in the z direction, the disk gravity de- 
creases while the projec tion of the radial halo gravity force 
in the z direction rises dSchulz fc Str uckl l200iri . This effect 
is particularly important in the case of low ram pressure 
regimes as the ones we are interested in here. The condition 



for stripping expressed in Equation (7) implicitly assumes 
instantaneous stripping from the disk and discrete ram pres- 
sure stripping events such that the loss of mass approaches 
an asymptotic limit. If instead ram pressure is a continuous 
process with time scales comparable with the orbital time 
of the satellite, the stripping radius calculated analytically 
underestimates the physical size of the gaseous disk. Such a 
continuous ram pressure stripping occurs here not because 
of turbulent or viscous stripping but because the gas at the 
edge of the disk heats up responding to the compression 
exerted by the outer medium and can therefore become un- 
bound by increasing its energy despite the ram pressure wind 
being constant. Timescales of both Kelvin-Helmoltz instabil- 
ities and artificial viscous forces are close to 10 Gyr because 
of, respectively, stabilization due to the cuspy halo potential 
of the LMC and the high mass resolution employed, and are 
therefore negligible over the timescales explored here. 

We estimate Rstr as the radius which contains 99% of 
the gas particles in the midplane of the disk. The mass of 
gas that remains in the disk can be computed summing the 
disk particles inside a cylinder of radius Rstr and thickness 
1 kpc. The mass of gas bound to the satellite at a time t is 
instead not directly related to Rstr, but is determined by the 
potential of the dark halo. We calculate it balancing kinetic, 
thermal and potential energy in concentric spherical shells 
centered on the LMC, defining the gas tidal radius Vg as the 
radius of the most distant bound shell. 

The mass of cold gas stripped from the LMC's disk nistr 
decreases as we decrease the mass of the halo particles: in 
HI nistr is ~ 1.6 times larger than in H3 (Table 2). Mass 
resolution also determinates the shape of the front edge and 
influences the morphology of the gaseous disk. Fig. 3 shows 
the gas density in the disk after 2 Gyrs for the HI, H2 and 
H3 runs. Numerical effects in the low resolution run are quite 
evident: the massive hot particles produce particle - parti- 
cle interactions with gas in the galaxy, creating holes and 
discontinuities along the front edge of the disk. The spiral 
pattern is partially destroyed and the disk looses its symme- 
try. Scattering reduces the stripping radius by a factor of 1.4 
and produces a sharp truncation, while the effect of pure ram 
pressure is much more gradual, since halo particles stream 
around the disk. Finally, the removed gas forms an asymmet- 
ric and turbulent tail. The high resolution H3 run appears 
very symmetric with a flat and uniform edge and a regu- 
lar stream, which flows around the disk and creates a low 
density regio n behind (Fig. 4). As previously noticed by sev- 
eral authors touilis. Moore fc Boweill2000l: lOuihs fc Moor'3 
2001; Schulz & Struck 2001), a bow shock forms in front of 
the disk, increasing the local ICM temperature from 10® to 
2.5 X 10® de grees. The intermediate run H2 is similar to HI, 
although the gas tail is more pronounced and the amount of 
stripped gas is less than 10% larger. The resolution (mass ra- 
tio of 1:2) of the H2 run is to be a good compromise between 
the necessity of solving the gas dynamics and the number of 
particles to use in a full production simulation. 

The runs LI, L2 and L3 have a gas density {ph = 
3.5 x 10~^cm~^) and satellite velocity (175 kms~^) simi- 
lar to our final LMC's orbit with apocentric distance ~ 50 
kpc and apo/peri ratio ~ 2.5 : 1 (see section 4.1). Reducing 
the ram pressure by a factor of three decreases the stripped 
mass by ~ 20%, while the amount of unbound gas is al- 
most 60% smaller. This indicates that a lower ram pressure 
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Table 2. LMC wind tunnel simulations. For each run the density pi^ of the hot medium and the velocity v 
of the satellite are indicated. The resolution of each model is expressed in terms of the mass ratio rrih/mdisk 
between halo and disk particles. N is the number of hot particles, while the last four columns give the ram 
pressure stripping radius Rstr and the mass stripped from the disk nigtr, the gas tidal radius Vg and the 
mass of gas unbound to the galaxy m^^i,. In the case of edge on galaxies the stripping radius -Rstr is the 
average value between the two sides of the disk. All the values are calculated at t = 2 Gyrs. 
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Figure 4. H3 run: density distribution of gas in a thin slice per- 
pendicular to the LMC disk. The satellite is moving face on to- 
ward the left side of the picture while a shock front forms in the 
external medium in front of the disk. 




Figure 3. From the top to the bottom: column density distribu- 
tion for the LMC test models HI, H2, H3. The scale is logarithmic 
with blue corresponding to the density of the hot halo gas. 



is still efficient in stripping gas from the disk, but part of the 
removed material is now bound to the satellite by the dark 
halo gravity. The simulations were computed using three dif- 
ferent resolutions, with mh/hdisk equal to 2, 1 and 0.5. The 
stripped mass converges for nih/mdisk < 1- 

The dependence of ram pressure on the satel- 
lite orientation is not a simple function of the an- 
gle between the d isk plane and the orbital velocity 
JVollmer et"al]|200ll) and can depe nd on the ram pressure 
^nd the structure of the satellite ("Abadi. Moore & Bowe,' 
1999; .Quilis, Moore fc Bowe r 2000; Sc hulz fc S truck 200 



i 



Marcolini. Brighenti fc D'Erc olc 2003). In the simulations 
El and E2 the galaxy is moving edge on through the hot 
gas. Comparing the amount of stripped mass of El and E2 
with the results of the corresponding face on runs H2 and 
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Figure 5. Orbital separation for the MW - LMC system 

LI, it is clear that the stripping is not very sensitive to 
disk orientation. This implies that the amount of gas in the 
Stream will only slightly depend on the initial inclination of 
the LMC disk and on the precession eventually induced by 
tidal forces from the MW. 



4 INTERACTING RUNS 
4.1 Orbital parameters 

The morphology of the Stream indicates that the orbit of the 
LMC is nearly polar and cou nter-clockwise, as seen from the 
Sun, while kinematical data (Ivan der Marel eral]l2002h im- 
ply that the present position at ~ 50 kpc from the Galactic 
Center, is close to the perigalacticon. Since the mass ratio 
between the MW and its satellite is about 50:1, we expect 
that the effects of the dynamical friction are not negligible 
on time scales comparable with the orbital period. On the 
other hand, at each pericentric passage the LMC loses mass 
from tidal stripping, increasing its decay time: dynamical 
friction becomes less effective as the satell ite approaches the 
Galaxy (fColpi. Maver fc Governatdll999ll . 

Dynamical friction implies that we can not pre- 
cisely make a backwards integration of the LMC's or- 
bit, therefore we make several low resolution test simula- 
tions (3 X 10^ particles) in order to find the correct start- 
ing point for the LMC 4 Gyrs in the past such that it 
ends up with the present distance, velocity and orienta- 
tion on the sky. We find orbital parameters close t o the 
ones obtained bv iGardiner. Sawa. fc Fuiimotol il994ft and 
ICardiner fc Noeuchil (Il996l) . In particular the orbital plane 
is perpendicular to the Galactic disk, the present eccentric- 
ity (defined as {rapo — rper)/{rapo + rper), with rapo and rper 
apogalacticon and perigalacticon distances) and the orbital 
period are respectively 0.44 and ~ 2 Gyrs. The orbital sep- 
aration between the MW and the LMC is plotted in Fig. 5 
for the past 4 Gyrs and 4 Gyrs into the future. The satel- 
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Figure 6. Final LMC rotation curve. 

lite is currently at 50 kpc from the Galactic center and is 
moving away from the perigalacticon (at 45 kpc) with a ve- 
locity of 250kms~^. From the plot it appears evident that 
the orbit is slowly sinking, with each apogalactic distance 
~ 20% smaller than the previous one. We will show in the 
next subsection that the change in the orbital parameters 
due to dynamical friction strongly affects the tidal forces 
and rate of ram pressure stripping. 

In choosing the initial inclination and position of the 
line of nodes we make the approximation that they do not 
change duri ng the interaction. In parti cular we adopt the 
values from Ivan der Marel et al.l i2002h : respectively i = 
34.7° and Q — 129.9°. This choice could partially influence 
the final stellar structure of the disk, but not the ram pres- 
sure stripping process, as seen in the previous section. 

4.2 The Stream 

The close interaction with the MW strongly perturbs the 
entire structure of the satellite. In particular during the 4 
Gyrs preceding the present epoch the LMC loses ~ 80% of 
its dark halo, and is currently tidally truncated at a radius 
of 11 kpc, consistent with the values provided by analytic 
calculations (JVoinbcrE 2000; van der Marel ct al. 2002) and 
with observations of the outer regions of the disk (|Ir-v^ril 
I1991I: Ivan der Marell 1200 ill . The contribution of the dark 
matter halo to the final rotation curve within 8 kpc from 
the center decreases by ~ 5 km s~^ and the disk component 
predominates in the inner region (Fig. 6). 

While all the stars within the tidal radius are likely to 
be gravitationally bound, this is not true for gas. In fact, 
part of the gas that lies within the tidal radius can have 
enough thermal energy to escape from the potential well of 
the galaxy. Fig. 7 illustrates the evolution of the baryonic 
mass as function of time. The dotted and solid curve refers 
to the steUar and gaseus component respectively while the 
dashed curve represents the gas within the tidal radius which 
is actually bound to the LMC. 
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Figure 7. Fraction of stars and gas bound to the satellite dur- 
ing the last 4 Gyrs. Both the gas within the tidal radius and 
the bound gas are plotted. Mi is the initial mass of each disk 
component. The vertical solid lines indicate the perigalactica. 




Figure 8. Present day distribution of stars from the LMC disk 
in a plane perpendicular to the Galactic plane. 



After two perigalactic passages the satellite looses 
roughly ~ 12% of its gas and ~ 4% of stars. A pure tidal 
stripping model would remove similar amounts of both stars 
and gas unless the initial gas disk was significantly more ex- 
tended. The larger amount of stripped gas is due primar- 
ily to ram-pressure stripping: the rate at which gas is lost 
increases as the LMC approaches perigalacticon (the solid 
vertical lines in the plot), corresponding to higher densities 
of the diffuse halo gas and to higher velocities of the satellite 
along the orbit. 

Stars are stripped from the LMC only during the last 
perigalacticon, since the previous perigalactic passage is too 
far from the MW to produce significative tidal shocks and 
to perturb the stellar structure of the satellite. The present 
stellar distribution (Fig. 8) indicates that a trailing and a 
small leading arm are forming (with surface density ~ 3 x 
lO~^M0pc~^), but most of the stars stripped from the disk 
are still bound to the LMC and form a large spheroidal 
component. 

Contrary to the stars, gas is stripped during the first 
passage at perigalacticon, forming a continuous stream that 
lies in a thin plane (width ~ 25 kpc) perpendicular to the 
disk of the MW (Fig. 9) and between 50 and 75 kpc from the 
Galactic center. Our simulations predict that the MS forms a 
great circle, but most of the gas is lost recently and it is now 
at an angular distance between 10° and 120° from the LMC 
(comparable with the range provided by observations). The 
mass of the stripped gas (1.4 x 10* M©) is of the same order 
of magnitude as the observed amount of neutral hydrogen 
in the Stream bv lPutman et al.l (jgOOSi) (2 x 10* Mq). 

Fig. 10 represents the projection of the Stream as seen 
from an observer situated at the Galactic center who is look- 
ing toward the South Galactic Pole. The logarithmic colour 
scale indicates the column density gradient which gradually 
decreases along the length of the Stream. At a fixed distance 
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Figure 9. Final distribution of gas from the LMC disk in a plane 
perpendicular to the Galactic plane 



from the LMC the structure of the Stream is not uniform but 
has clumpy and high density regions. The decrease in den- 
sity is by nearly two orders of magnitude along the Stream 
(Fig. 11), due to the combination of two factors. The gas 
close to the head of the Stream is pulled from the disk dur- 
ing the last perigalacticon, when the stripping rate is maxi- 
mum, whereas the tip is ~ 1 Gyr old and corresponds to an 
apogalactic passage. Moreover, the tail is formed by the first 
material stripped from the outer and low density regions of 
the exponential disk. This decrease in column density is a 
remarkable success for the ram-pressure scenario since tidal 
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Figure 10. Hydrogen column density along the simulated stream. 
The intensity values are on a logarithmic scale with white corre- 
sponding to the highest density regions within the LMC disk and 
dark blue to njj ~ 10^*cm~^. 



models generically produce streams with surface densities 
that fall off much more slowly. 

The difference in the characteristic time scales of hydro- 
dynamical and gravitational forces explains why a gaseous 
leading arm does not form through tidal stripping: gas is 
removed quite early from the outer regions of the LMC disk 
that during the second passage at the perigalacticon starts 
forming a leading arm feature. The great circle of stripped 
gas actually appears on the sky as a leading stream and 
could in principle reproduce observations of the leading arm 
(Fig. 9). This material was stripped from the satellite dur- 
ing the first orbit and is falling to the Galactic center in 



10=1 r , , , , 1 , , , , \ , , q 




distance (degree) 

Figure 11. Hydrogen column density along the Stream. The sys- 
tem of coordinates is chosen in order to have the LMC at de- 
grees. 

the Northern Galactic hemisphere, intersecting the actual 
position of the satellite. 

The total amount of unbound gas after 4 Gyrs is com- 
parable with the gas stripped in one orbital time from an 
LMC model in a test tube simulation with densities and ve- 
locities characteristic of a perigalactic passage. This means 
that the tidal forces acting on the satellite contribute to ram 
pressure processes: as a consequence of the resulting asym- 
metrical potential and distorted stellar disk, it is easier to 
strip more gas from the LMC's disk since the gravitational 
restoring force is weaker. Thus even a low density gaseus 
Galactic halo is able to remove a significant amount of gas 
from the LMC. 

4-2.1 Adiabatic versus cooling runs 

We have presented results from simulations that do not in- 
clude radiative cooling. As discussed in section 3, in this 
case the gas in the disk of the LMC will increase its tem- 
perature as it is strongly compressed by the outer medium. 
In the adiabatic regime the gas can only lose energy by ex- 
pansion, which requires a long timescale, and can become 
easily unbound due to this input of thermal energy. We have 
run a simulation with radiative cooling adopting a standard 
cooling function for primordial gas (hydrogen and helium). 
Since the gas quickly heats up to ~ 10^ K, the peak of the 
cooling curve, then cools down on a timescale much shorter 
than the dynamical time, shrinking in radius compared to 
its original location. As a result a smaller gaseous disk re- 
sults after a few cycles of heating and rapid cooling which 
sits deeper in the potential well and as a result no gas is 
stripped. However, what we are witnessing is probably the 
result of overcooling since we are neglecting the effect of 
star formation, namely heating from supernovae feedback 
and radiation from massive stellar associations. The LMC 
in particular shows widespread intense star formation across 
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its disk and bar. In reality a multi-phase ISM will result that 
will be eventually close to the condition of pressure equilib- 
rium with the outer medium and will not be able to cool so 
efficiently. Until we can effectively model this multi-phase 
medium, we believe that the results of the adiabatic simu- 
lations are closer to reality. 



4.3 The LMC disk 

The Galactic disk remains unperturbed during the inter- 
action, while the morphology of the satellite changes drasti- 
cally. In particular we will focus on the changes in the stellar 
and gaseous structure of the LMC. 

After 4 Gyrs from the beginning of the simulations al- 
most the entire stellar component is still bound to the satel- 
lite, but less than 90% of the stars lies in the thin disk, which 
is surrounded by a stellar spheroid with a radius of ~ 11 
kpc (Fig. 12). The structure of this spheroidal component is 
quite complex, since a strong warp that wraps 180° around 
the LMC is superimposed on the underlying low density 
stellar distribution. The warp forms after the first passage 
at the perigalacticon and involves a large fraction of the stel- 
lar particles. It extends up ~ 5 kpc out of the thin disk and 
is aligned with the ma jor axis of the sat ellite. The existence 
of a warp is inferred bv lKim et alJ il998[l from the gas veloc- 
ity fie ld in the disk of the LMC, bv Ivan der Marel fc Cionil 
J200ll) and lOlsen fc Salvkl ll2002f) using the apparent mag- 
nitude of different classes of stars as a distance indicator. 
The diffuse spheroidal component has an average surface 
density of ~ 5 x lO~^M0pc~^, which becomes one order 
of magnitude larger if we consider only the warp region. 
The warp deformation is partially directed along the line 
of sight and could have implications for the microlensing 
results, increasing the optical dep th caused by self-lensing 
JZaritskv fc Linll997l : lzhao fc Ev ans 2000). The LMC's disk 
appears elongated, with a mean ellipticity e — 0.38 (where 
e = 1 — q and q is the axial ratio), in a direction roughly 
aligned with the Galactic center and perpendicular to the 
Stream. Th ese values agree very well with the observations 
provided bv lvan der Marel i200l|) . Within the inner 6-7 kpc 
of the model LMC a strong and asymmetric bar forms after 
the first perigalacticon. 

The isodensity contours of the face-on disk were fitted 
with ellipses using the MIDAS task FIT/ELL3. The eUipse 
centers, ellipticity and major axis position angle were con- 
sidered as free parameters of the fitting procedure. The re- 
sults of this projected analysis are illustrated in Fig. 13. The 
first plot on the top (solid line) shows the B-band surface 
brightness profile of the present face-on LMC, for which the 
best fitting exponential curve has a scale length of 1.5 kpc. 
The dashed line is the original exponential disk, with a scale 
length of ~ 1.7 kpc. The second and the third plot illustrate 
the variation of the position angle and eccentricity as a func- 
tion of radius. The position angle is almost constant in the 
central 4 kpc of the disk and presents a sharp twist beyond 
this radius, which corresponds to the transition between the 
bar region and the outer density contours characterized by 
ellipticities <0.2. At radii > 8 kpc the disk appears more 
elongated, with a further increase in ellipticities and a twist 
of ~ 50° in the position angle. The last plot on the bottom 
indicates the drift of the countour centers in the plane of the 




Figure 12. LMC stellar surface density. Colours represent a loga- 
rithmic scale where white corresponds to a density of lO^M0pc~^ 
and blue to 10~^Mqpc~^. From the top to the bottom: initial 
conditions, final face on and edge on projection. The bar and the 
strong warp are visible. 
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Figure 13. Photometric analysis of a face-on projection of the 
LMC disk. From the top to the bottom: B-band surface bright- 
ness, major axis position angle and ellipticity profiles. The last 
plot on the bottom represents the drift of the ellipse center in the 
plane of the disk. 
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Figure 14. Projected velocity (left) and velocity dispersion 
(right) profiles for three different orientations of the LMC disk. 
From the top to the bottom: face-on, edge-on and end-on. 



disk, which moves consistent with observations by ahnost 1 
kpc. 

Fig. 14 shows the kinematic profiles of three different 
orientations of the LMC disk observed today: from the top 
to the bottom face-on, edge-on and end-on projections. The 
peak in the velocity dispersion in the central 6-8 kpc is 
mainly due to the bar formation, while the external regions 
appear to be dominated by an hot tidally perturbed com- 
ponent with (Tp ~ 20 km s^'^. 

The gaseous disk is completely dominated by stripping 
processes (Fig 15). It retains a symmetric structure only 
within the central 12 kpc, in connection with the thin com- 
ponent of the stellar disk. Within this radius we do not ob- 
serve a significative displacement between the center of the 
stellar distribution and the gaseus disk. The ram pressure 
stripping radius after 4 Gyrs is ^ 6 kpc, a factor of three 
smaller than the initial disk radi us and compatible with the 
results provided by HIPASS data jPutman. Gibson fc SmithI 
ITqQSi). 



5 CONCLUSIONS 

We have carried out high resolution gravita- 
tional/hydrodynamical simulations of the interaction 
between the LMC and Milky Way. As the LMC spirals 
inwards towards the Galaxy it suffers ever increasing 
gravitational tidal forces and hydrodynamical stripping. 
Our simulations cover the previous 4 Gyrs of the orbit of 
the LMC such that at the final time it ends up at its correct 
location within the Galactic halo and with the observed 
inclination. We find that the combined effects of gravity 
and ram-pressure stripping can account for the majority of 
the LMC's kinematical and morphological features and the 
morphology of Magellanic Stream. 

Our model of the hot gaseous halo places a total of 
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Figure 15. LMC gaseous column density. Colours represent a log- 
arithmic scale where white corresponds to a density of lO'^^cm"'^ 
and violet to lO^^cm"^. From the top to the bottom: initial con- 
ditions, final face on and edge on projection. The bar and the first 
part of the Stream are visible. 



IQ Mq of ionised hydrogen within the dark matter halo of 
the Galaxy. The density at 50 kpc is « 8 x 10~^ atoms/cm^. 
This is sufficient to remove over 1O*M0 of gas from the 
LMC, close to the observed mass of the Stream. The stripped 
gas forms a great circle, or a polar ring, around the Galaxy, 
consistently with the recently discovered extension of the 
Stream into the Northern Hemisphere bv lBraun fc Thilkerl 
12004*). Less material was stripped during the pericentric 
passage 4 Gyrs ago since dynamical friction has moved the 
Clouds closer today. Thus we can reproduce the observed 
decrease in column density along the Stream. Very few stars 
are tidally removed from the LMC, but its disk becomes 
severely warped due to the tidal interaction. This creates 
a diffuse halo of stars within the LMC that may have been 
observed as self-lensing events. The disk of stars is elongated 
by 2:1 oriented towards the Galactic centre in agreement 
with the observations of van dcr Marcl (2001). 

Several improvements to our simulations could be made 
in order to study this fascinating interacting system in more 
detail. More detailed modelling of the star-formation and 
internal ISM of the LMC is necessary in order to study the 
role of gas stripping on star-formation, and the role of star- 
formation and the multiphase ISM on the stripping process. 
Including the interaction with the SMC which could sup- 
ply a significant amount of gas to the Stream and may be 
responsible for the Magellanic bridge and the leading arm. 
Including a more realistic treatment of the gas physics in or- 
der to study the detailed structure of the stripped gas and 
its interaction with the halo gas as it infalls through the 
galactic halo. 
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